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Background: 

Over  the  past  few  years  the  Naval  Research  Laboratory  has  investigated  the  use  of 
the  AN/SQS-53B  ASW  sensor  as  a  research  tool  in  seafloor  characterization  and 
mapping.  This  program  has  resulted  in  the  development  of  a  data  acquisition  system 
which  can  be  attached  to  the  sonar  receiver  to  obtain  digital  recordings  of  the  pulse 
returns.  The  latest  version  of  the  system  was  installed  on  the  USS  Monterey  (CG-61) 
prior  to  its  planned  transit  of  the  Red  Sea  in  November  1995. 

The  AN/SQS-53B  has  a  number  of  operating  modes  designed  to  optimize  its 
primary  mission  of  detecting  submarines  in  the  deep  ocean.  The  modes  of  most  interest 
in  acoustic  seafloor  mapping  are  the  Convergence  Zone  (CZ),  Bottom  Bounce  (BB),  and 
Bottom  Bounce  /  Tracking  (BB/TK)  modes.  In  the  first  two,  the  received  reflected 
acoustic  pulse  is  formed  into  twelve  nominally  conical  beams  each  of  which  is 
approximately  10°  wide  in  the  horizontal  plane.  The  receiver  beams  can  be  steered  so 
that  the  entire  120°  coverage  is  centered  on  any  azimuth.  The  beams  are  stabilized  so  that 
once  the  sector  center  of  the  120°  coverage  is  set,  each  beam  remains  pointed  in  a 
constant  direction,  relative  to  true  north,  regardless  of  the  ship’s  course.  The  beams  can 
also  be  steered  vertically  down  to  45°  depression  angle.  The  BB/TK  mode  is  similar, 
with  the  exception  that  the  projector  beam  is  constrained  to  a  30°  horizontal  coverage 
directed  along  the  sector  center.  This  characteristic  provides  increased  power  in  the 
direction  of  a  previously  located  target  for  tracking. 

The  horizontal  stabilization  of  the  sonar  beams  (built  into  the  on  board  electronic 
beamforming  and  the  ship's  inertial  navigation)  provides  a  capability  for  acoustic  image 
mapping  not  found  on  conventional  side-scan  sonars  or  imagery  derived  from  multibeam 
systems.  In  conventional  side-scan  systems  the  beam  look-direction  is  always 
perpendicular  to  the  tow  ship  track.  The  availability  of  multiple  beams,  each  stabilized  in 
a  different  geographic  look-direction,  allow  the  construction  of  an  acoustic  backscatter 
image  for  each  look  direction.  In  many  areas  of  the  world,  backscatter  from  the  seafloor 
is  highly  dependent  on  the  azimuthal  aspect  angle  of  the  incident  acoustic  signal.  This 
fact  has  historically  hindered  the  full  exploitation  of  conventional  side-scan  sonar  systems 
in  bottom  mapping.  An  example  of  the  aspect  dependence  of  the  backscatter  is  in  areas 
where  the  seafloor  is  highly  delineated  (such  as  in  the  Red  Sea  or  the  Mid-Atlantic 
Ridge).  When  the  linear  seafloor  features  are  illuminated  broadside,  the  acoustic  return  is 
much  higher  than  when  the  illumination  is  along  the  linear  direction.  The  development  of 
an  acoustic  imagery  mapping  product  requires  a  means  of  accounting  for  these  large 
aspect  angle  dependencies.  One  solution,  suggested  here,  is  that  separate  maps  be 
prepared  for  each  available  look  direction.  Such  a  data  set  would  then  prove  useful  in 
developing  methods  of  combining  look  directions  to  categorize  the  backscatter. 

The  Naval  Oceanographic  Office  conducted  a  multibeam  bathymetry  survey  in 
the  Red  Sea  during  the  same  time  the  Monterey  made  its  transit.  The  combination  of  this 
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data  with  the  AN/SQS-53  imagery  data  allows  increased  exploitation  of  both  data  sets. 
This  fusion  is  possible  because  of  the  precise  location  of  the  AN/SQS-53  based  on  GPS. 

Figures  1  and  2  show  the  location  of  an  area  of  the  Red  Sea  (Area  3  on  the  map) 
chosen  for  a  preliminary  analysis.  The  background  bathymetry  is  at  a  2-arc  minute 
resolution  (DBDBV)  and  the  Naval  Oceanographic  Office  multibeam  data  was  process  to 
a  100  meter  resolution  as  was  the  sonar  imagery. 

Approach: 

The  NRL  data  acquisition  system  digitized  data  from  each  of  the  twelve  receiver 
beams  formed  in  the  three  modes.  For  each  beam,  data  was  acquired  from  two  separate 
points  in  the  receiver.  These  points  were  essentially  before  and  after  the  Automatic  Gain 
Control  (AGC)  amplifier  but  after  the  Time  Varied  Gain  (TVG)  amplifier  in  the  receiver 
circuit.  These  points  are  designated  Pre-AGC  and  Post-AGC  channels.  Thus,  24 
channels  of  data  were  recorded.  Previous  NRL  experiments  indicated  that  Pre-AGC  data 
was  preferable  to  the  Post-AGC  data  for  bottom  mapping. 

Two  sonar  acoustic  data  sets,  Pre-  and  Post-AGC,  were  digitized  at  a  rate  of  250 
Hz  per  chaimel  giving  a  range  sampling  distance  of  3  meters.  The  data  sets  were  recorded 
on  optical  disks  in  a  "band-interleaved-by-pixel"  format.  Each  ping  or  scan  was  preceded 
by  a  header  containing  ping  number  and  time.  Other  fields  in  the  ping  header  were  used 
for  additional  information  including  navigation  and  sonar  operating  modes.  During  this 
deployment,  however,  the  latter  data  was  not  readily  available  for  digital  recording.  As  a 
result,  navigation  was  recorded  separately  (by  hand,  every  5  minutes)  using  a  GPS 
receiver.  Time  differences  between  the  GPS  and  the  computer  controlling  the  data 
acquisition  were  noted.  In  addition,  sonar  modes  and  sector  center  locations  were 
recorded  in  the  log  book. 

In  order  to  produce  an  acoustic  image  map  from  the  acquired  data,  it  is  necessary 
to  locate  the  position  of  the  ship  at  the  time  the  projector  transmit  ping  is  made.  This 
position,  combined  with  the  known  direction  of  the  beam  and  the  water  depth,  allows  the 
received  ping  to  be  geo-rectified  on  a  map  datum.  Subsequent  pings  can  thus  be 
combined  into  a  geo-rectified  image.  Due  to  the  ping  length  and  beamwidth 
combinations,  the  "raw"  imagery  has  pixels  with  large  differences  in  the  extent  of 
transmit  direction  and  cross  beam  direction.  A  separate  image  can  be  formed  fi-om  each 
beam.  As  noted,  the  image  is  analogous  to  a  side-scan  sonar  image  except  that  the  look 
directions  are  constant  (relative  to  true  north)  and  the  beams  are  not  all  directed  at  90 
degrees  to  the  ship's  course. 

The  initial  step  in  the  data  analysis  consists  of  de-multiplexing  the  time  sequence 
for  each  beam  and  merging  with  the  navigation  and  log  data .  This  was  accomplished 
based  on  the  time  of  the  ping  and  the  time  of  the  GPS  fix.  The  latter  was  accurate  to  1 
microsecond. 
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Twelve  separate  "channels,"  or  data  files,  were  constructed  for  each  ping.  Each 
file  contained  a  proto-image  that  was  swept  out  as  the  ship  moved  along  track.  The 
approach  to  geo-rectification  of  each  file  is  based  on  the  geometry  of  a  geo-stabilized  and 
a  relatively  constant  ship  course.  This  approach  has  proved  imsatisfactory  for  the  beams 
pointed  in  the  direction  of  the  ship's  heading  but  is  well  suited  to  the  other  beams.  The 
details  of  the  geo-rectification  and  the  "slant  range"  correction  using  a  flat  bottom 
assumption  are  included  in  the  Appendices.  The  procedure  consists  of  placing  each  ping 
(considered  as  a  line)  in  its  proper  geographic  position  based  on  the  interpolated  position 
of  the  start  of  the  ping  and  the  number  of  samples  to  the  end  giving  the  range  to  the  end 
of  the  ping  (assuming  a  constant  sound  velocity).  The  pings  are  thus  aligned  in  a 
rectangular  box  by  padding  out  the  beginnings  and  ends  with  zeros.  The  location  of  the 
comers  of  the  boxes  are  known.  The  aligned  ping  data  is  then  resampled  in  a  direction 
perpendicular  to  the  ping  direction,  resulting  in  a  file  with  equally  spaced  data  in  both 
directions.  The  "gridded"  data  is  imported  into  an  image  processing  computer  program 
(ERDAS  Imagine  ver  8.1)  and  the  geo-rectification  completed. 

Once  the  data  is  geo-rectified  and  put  in  image  format,  it  can  be  compared  to 
bathymetry  and  other  data;  a  number  of  image  enhancements  can  be  performed  to  aid  in 
data  interpretation. 

Results: 

The  processing  steps  described  above  are  shown  in  Figures  3  through  5c.  Figure 

3  shows  the  raw  data  for  a  single  beam.  This  particular  data  set  contains  pings  #505 
through  #600  of  file  27  (a  Post-AGC  set).  Figure  4  shows  the  "gridded"  data  for  all 
twelve  beams.  Figures  5a  and  5b  (different  scales)  show  the  geo-rectified  Post-AGC 
imagery  and  Pre-AGC  imagery,  respectively,  for  a  look  direction  of  52°.  In  Figures  5a- 
5c,  the  ship’s  track  is  from  lower  right  to  upper  left.  The  ping  origins  are  along  the  left 
(lower)  edge  of  the  data.  The  bright  bands  on  the  right  (upper)  edge  of  the  image  are 
artifacts  that  have  not  yet  been  corrected.  Comparison  of  Figures  5a  and  5b  reveal 
differences  in  texture,  i.e.,  some  bottom  features  show  up  better  in  one  or  the  other  of 
these  images.  For  example,  the  highlight  returns  in  the  lower  right  of  Figure  5b  originate 
from  a  seamount;  however,  in  Figure  5a  only  the  shadow  of  the  seamount  is  evident  (the 
highlights  are  obscured  by  the  overall  bottom  reverberation).  Figure  5c  shows  that  the 
best  aspects  of  both  the  Post-AGC  (red)  and  Pre-AGC  (green)  images  are  enhanced  when 
the  images  are  merged. 

It  is  difficult  to  display  all  twelve  geo-rectified  beams  because  they  overlap  each 
other.  Figure  6a  shows  only  ten  look  directions  superimposed,  but  illustrates  the 
difficulty  of  “seeing”  multiple  aspect-angles.  For  comparison.  Figure  6b  shows  the 
details  of  one  look  direction  at  52°.  However,  on  the  computer  it  is  possible  to  flicker 
between  two  or  more  beams  by  varying  the  transparency.  One  method  of  simultaneously 
viewing  up  to  three  beams  is  to  use  color  coding.  In  Figure  7a,  three  beams  with  look 
directions  of  2  °  (upper),  52  °  (center)  and  92  °  (lower)  are  superimposed.  Inside  the 
boxed  area  (Figure  7a),  the  3  beams  have  been  converted  to  color  (Figure  7b)  by 
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assigning  a  different  color  gun  to  each  beam:  red  for  92°,  green  for  52  °,  blue  for  2  °.  For 
features  common  to  all  three  beams  the  resultant  color  should  be  white  (if  the  magnitudes 
are  equal).  If  a  feature  occurs  in  only  one  beam  then  it  will  appear  in  the  color  assigned 
to  that  beam.  As  can  be  seen  in  Figure  7b,  many  of  the  ridge  features  are  not  common  in 
all  the  beams.  This  is  attributed  to  the  different  aspect  angles.  The  use  of  color  in  this 
type  of  analysis  requires  further  investigation. 

Comparing  the  separate  beams  with  bathymetry  can  be  difficult.  One  method  is 
to  drape  the  imagery  over  the  bathymetry  and  display  the  result  in  perspective  shaded 
relief  A  segment  of  the  bottom  bathymetry  in  the  vicinity  of  image  ni27/700/52  (see 
Figure  1  for  location)  is  shown  in  Figure  8  in  shaded  relief  Figure  9  shows  the  result  of 
draping  the  imagery  (m27/700/52)  on  the  bathymetry.  (In  Figures  8  and  9  we  are  looking 
north).  This  overlay  demonstrates  that  while  there  is  a  strong  correlation  of  the 
backscatter  with  many  bottom  features,  this  is  not  always  the  case. 

The  sonar  imagery  in  this  area  has  also  been  compared  to  previously  collected 
GLORIA  data.  GLORIA  is  a  3.5  kHz  side-scan  sonar  that  has  been  used  world  wide  in 
acoustic  bottom  imaging.  Figure  10  is  a  portion  of  the  GLORIA  data  (a  photo  copy  from 
a  journal  article)  that  has  been  geo-rectified  after  scanning.  Preliminary  analysis  indicates 
that  although  the  GLORIA  along-track  resolution  is  much  better  than  that  of  the 
AN/SQS-53B,  the  AN/SQS-53B  detected  all  the  linear  features  that  are  evident  in  the 
GLORIA  data. 

The  performance  of  the  AN/SQS-53B  in  shallow  water  (less  than  100  meter)  was 
also  investigated  using  the  available  data  set.  This  data  is  shown  in  Figures  1 1  and  12a 
and  12b  for  the  Strait  of  Bab  al  Mandab  that  connects  the  Red  Sea  and  the  Gulf  of  Aden. 
The  sonar  signals  in  Figures  1 1  and  12  (a  and  b)  use  slightly  different  beam 
combinations.  Figures  1 1  shows  beams  at  angles  of 270  °  and  10  °,  and  Figures  12a  and 
12b  shows  beams  at  angle  of  10  °.  Figure  12a  shows  the  acoustic  returns  from  the 
shoaling  areas  around  Perim  Island.  A  bathymetric  chart  of  this  area  was  digitized  and 
the  acoustic  imagery  of  Figure  12a  superimposed  on  the  chart  (Fig.  12b).  In  Figure  12b 
strong  returns  are  shown  in  black,  which  is  the  reverse  of  Figure  12a  (and  all  the  other 
figures).  This  is  done  to  allow  low  returns  to  be  “transparent .”  This  overlay  is  much 
easier  to  interpret  on  the  computer  display  since  the  transparency  of  the  imagery  can  be 
readily  controlled.  A  detailed  comparison  of  the  chart  and  the  image  shows  a  close 
correlation  with  all  the  charted  features.  However,  some  returns  on  the  sonar  imagery  are 
not  charted.  This  demonstrates  the  value  for  detecting  possible  hazards  to  navigation  and 
for  delineating  bottom  structure  at  high  spatial  resolution  (compared  to  the  chart)  in 
shallow  water.  The  data  was  obtained  as  the  ship  transited  the  Strait  of  Bab  al  Mandab  at 
a  speed  of  25  knots.  The  capability  demonstrated  here  suggests  the  use  of  this  sonar  and 
geo-referencing  for  aiding  the  planning  and  execution  of  detailed  and  accurate 
hydrographic  surveying. 

The  particular  geo-rectification  technique  used  up  to  this  point  does  not  show 
returns  from  directly  in  front  of  the  ship  well.  This  is  illustrated  in  Figure  13  which 
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shows  the  beam  with  a  look  direction  of  320”.  In  the  lower  segment  a  fairly  good  image 
results  from  the  ship  track  being  about  340”.  However,  when  the  ship  turns  left  to  a 
course  of  about  3 1 5"  the  resultant  image  is  poor.  The  fix  for  this  situation  is  a  "full-up" 
geo-rectification  for  each  pixel  on  the  bottom  where  the  cumulative  effects  of  all  pings 
incident  from  a  particular  aspect  angle  are  accounted  for  and  the  beam  pattern  of  the 
sonar  is  included  in  the  calculation.  This  is  a  more  computer  intensive  method  than  the 
one  used,  but  can  be  expected  to  give  superior  results.  Alternate  means  of  handling  multi¬ 
aspect  acoustic  imagery  are  being  examined,  including  3-color  composites  and  controlled 
transparency.  Further  efforts  in  this  area  are  needed  to  fully  exploit  the  information  in 
the  imagery. 

Summary: 

Preparation  of  geo-rectified  AN/SQS-53B  acoustic  imagery,  using  data  collected 
during  November  and  December  1995  transits  of  the  Red  Sea,  has  been  completed  and 
superimposed  with  Naval  Oceanographic  Office  multibeam  bathymetry  collect  during  the 
same  period.  Preliminary  comparison  of  AN/SQS  53B  sonar  data  with  the  multibeam 
bathymetry  shows  positive  correlations  between  backscatter  and  seafloor  relief  (i.e.,  slope 
aspect)  in  both  deep  and  shallow  water.  High  reverberation  levels  are  associated  with 
deep-water  linear  ridges  when  insonified  normal  to  their  strike.  In  contrast,  low  levels 
result  when  the  look  direction  parallels  the  lineations.  In  addition,  features  associated 
with  island  shoals  were  evident  in  reverberation  data  collected  at  a  distance  of  4  miles 
from  the  island.  It  is  demonstrated  that  the  AN/SQS-53B  acoustic  backscatter  data  can  be 
collected,  geo-rectified  and  archived  into  a  potentially  valuable  product. 

Conclusions: 

The  AN/SQS-53B  is  the  primary  active  sonar  submarine  detection  equipment  in 
the  Navy  for  surface  combatants  and  is  installed  in  numerous  ships  including  cruisers, 
destroyers  and  frigates.  The  normal  deployment  of  these  ships  could  allow  the  collection 
of  a  huge  amount  of  acoustic  seafloor  backscatter  data  over  time.  A  digital  acoustic 
imagery  product,  populated  by  operational  AN/SQS-53B  sonars,  can  be  developed  to 
fully  exploit  the  capabilities  of  this  sonar.  In  addition,  on  board  processing  can  be 
developed  to  allow  the  direct  generation  of  geo-rectified  imagery  for  comparison  to 
archived  NIMA  data.  The  advent  of  accurate  GPS  positioning  for  the  ship  combined  with 
the  computing  power  required  for  image  processing  has  made  the  concepts  outlined  here 
feasible. 
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Figure  1.  AN/SQS-53B  sonar  imagery  (black  and  white)  superimposed  onto  multibeam  bathymetry 
(Color).  Bottom  depths  range  from  shallow  water  (red)  to  deep  water  (dark  blue).  Solid  circles  mark  the 
northward  course  of  the  ship  and  open  circles  mark  the  southward  course. 

Figure  2.  Enlargement  of  Area  3  shown  in  Figure  1.  Bottom  depths  range  from  shallow  water  (red)  to  deep 
water  (dark  blue).  Arrows  indicate  the  ship’s  course  and  heading,  and  origin  of  the  outgoing  pings. 

Figure  3.  Single-ping,  AN/SQS-53B  raw  sonar  data  for  ping  numbers  505-606. 

Figure  4.  Gridded  AN/SQS-53B  sonar  data  for  all  twelve  beams. 

Figure  5a.  Geo-rectified,  Post-AGC,  AN/SQS-53B  imagery  for  052"  look  direction.  The  ship’s  heading 
(relative  to  the  box)  is  from  lower  right  to  upper  left  and  the  ship’s  track  coincides  with  the  left  (lower)  side 
of  the  imagery.  For  reference  with  Figure  2,  the  coordinate  for  the  lower  right  comer  of  the  box  is 
17°N/41°R 

Figure  5b.  Geo-rectified,  Pre-AGC,  AN/SQS-53B  imagery  for  052°  look  direction.  The  ship’s  heading  is 
from  lower  right  to  upper  left  and  the  ship’s  track  coincides  with  the  left  (lower)  side  of  the  imagery. 

Figure  5c.  Geo-rectified,  merged  Pre-  and  Post-AGC,  AN/SQS-53B  imagery  for  052P  look  direction. 

Figure  6a.  Geo-rectified,  AN/SQS-53B  sonar  imagery  for  ten  look  directions.  The  ship’s  heading  is  from 
lower  right  to  upper  left  and  the  ship’s  track  coincides  with  the  left  (lower)  side  of  the  imagery. 

Figure  6b.  Geo-rectified  AN/SQS-53B  sonar  imagery  for  approximately  the  same  field  shown  in  Figure  6a, 
but  for  a  look  direction  of  primarily  052°  (rectangular  area).  Portions  of  the  002°  and  092°  look  directions 
show  at  the  top  and  bottom  edges  of  the  rectangle,  respectively.  The  ship’s  heading  is  from  lower  right  to 
upper  left  and  the  ship’s  track  coincides  with  the  left  (lower)  side  of  the  imagery. 

Figure  7a.  Geo-rectified  AN/SQS-53B  sonar  imagery  for  002°  (upper),  052°  (center)  and  092°  (lower)  look 
directions  (field  is  essentially  the  same  as  in  Figures  6a  and  6b).  Rectangular  area  shows  a  composite  of  the 
three  beams. 

Figure  7b.  Enlargement  of  rectangular  area  in  Figure  7a.  The  colors  represent  the  different  geo-rectified 
beams  or  look  directions:  red  for  092°,  green  for  052°  and  blue  for  002°. 

Figure  8,  Three-dimensional,  shaded  bathymetry  in  the  vicinity  for  image  M27/700/52  (looking  north). 

Figure  9.  AN/SQS-53B  sonar  image  M27/700/52  draped  over  shaded  bathymetry  shown  in  Figure  8 
(looking  north). 

Figure  10.  Geo-rectified  GLORIA  data  (3.5  kHz)  for  the  same  area  shown  in  Figure  5a. 

Figure  11.  AN/SQS-53B  shallow  water  performance  at  beam  angles  of  OKF  and  270°.  Data  collected  at  25 
kts. 

Figure  12a.  AN/SQS-53B  shallow  water  performance  at  beam  angle  of  ICP  showing  acoustic  returns 
(bright  areas)  from  the  shoal  areas  surrounding  Perim  Island. 

Figure  12b.  Acoustic  returns  shown  in  Figure  12a  (but  shown  here  as  dark  areas)  superimposed  on 
bathymetric  chart  of  the  area  surrounding  Perim  Island. 
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Figure  13,  Geo-rectified  AN/SQS-53B  sonar  image  forward  of  the  ship  as  it  heads  in  a  northward  direction. 
The  track  of  the  ship  is  traced  by  the  bright  or  whited-out  portion  of  the  imagery.  Note  that  when  the  ship 
changed  course,  the  look  direction  shifted  from  the  port  side  to  the  starboard  side  of  the  vessel  due  to  the 
receiver  beams  being  stabilized  (i.e,,  each  beam  remains  pointed  in  a  constant  direction,  relative  to  true 
north,  regardless  of  the  ship’s  course). 
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Appendices 


1 .  Conversion  from  range  to  horizontal  distance  under  a  flat  bottom  assumption. 

2.  Specialized  AN/SQS-53B  Geo-rectiflcation. 
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AN/SQS-53B  /  Red  Sea 


Bathymetry  and  Sonar  coverage 


Figure  1.  AN/SQS-53B  sonar  imagery  (black  and  white)  superimposed  onto  multibeam  bathymetry  (Color). 
Bottom  depths  range  from  shallow  water  (red)  to  deep  water  (dark  blue).  Solid  circles  mark  the  northward 
course  of  the  ship  and  open  circles  mark  the  southward  course. 


AN/SQS-53B  Sonar  / 
Red  Sea  Bathymetry 
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Figure  2.  Enlargement  of  Area  3  shown  in  Figure  1 .  Bottom  depths  range  from  shallow  water  (red)  to  deep 
water  (dark  blue).  Arrows  indicate  the  ship’s  course  and  heading,  and  origin  of  the  outgoing  pings. 


Figure  3.  Single-ping,  AN/SQS-53B  raw  sonar  data  for  ping  numbers  505-606. 
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Figure  4.  Gridded  AN/SQS-53B  sonar  data  for  all  twelve  beams. 


Figure  5a.  Geo-rectified,  Post-AGC,  AN/SQS-53B  imagery  for  052°  look  direction.  The  ship’s  heading 
(relative  to  the  box)  is  from  lower  right  to  upper  left  and  the  ship’s  track  coincides  with  the  left  (lower)  side  of 
the  imagery.  For  reference  with  Figure  2,  the  coordinate  for  the  lower  right  comer  of  the  box  is  17°N/41°E. 


fled,  Pre-AGC,  AN/SQS-53B  imagery  for 
eft  and  the  ship’s  track  coincides  with  the 


Figure  5c.  Geo-rectified,  merged  Pre-  and  Post-AGC,  AN/SQS-53B  imagery  for  052°  look  direction. 


Reset  Window  Tools 


Figure  6a.  Geo-rectified,  AN/SQS-53B  sonar  imagery  for  ten  look  directions.  The  ship’s  heading  is  from  lower 
right  to  upper  left  and  the  ship’s  track  coincides  with  the  left  (lower)  side  of  the  imagery. 


Figure  6b.  Geo-rectified  AN/SQS-53B  sonar  imagery  for  approximately  the  same  field  shown  in  Figure  6a,  but 
for  a  look  direction  of  primarily  052°  (rectangular  area).  Portions  of  the  002°  and  092°  look  directions  show  at 
the  top  and  bottom  edges  of  the  rectangle,  respectively.  The  ship’s  heading  is  from  lower  right  to  upper  left  and 
the  ship’s  track  coincides  with  the  left  (lower)  side  of  the  imagery. 
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Figure  7a.  Geo-rectified  AN/SQS-53B  sonar  imagery  for  002°  (upper),  052°  (center)  and  092°  (lower)  look 
directions  (field  is  essentially  the  same  as  in  Figures  6a  and  6b).  Rectangular  area  shows  a  composite  of  the 
three  beams. 


Figure  7b.  Enlargement  of  rectangular  area  in  Figure  7a.  The  colors  represent  the  different  geo-rectified  beams 
or  look  directions:  red  for  092°,  green  for  052°  and  blue  for  002°. 


Figure  8.  Three-dimensional,  shaded  bathymetry  in  the  vicinity  for  image  M27/700/52  (looking  north). 
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Figure  9.  AN/SQS-53B  sonar  image  M27/700/52  draped  over  shaded  bathymetry  shown  in  Figure  8  (looking 
north). 


Figure  10.  Geo-rectified  GLORIA  data  (3.5  kHz)  for  the  same  area  shown  in  Figure  5a. 


AN/SQS-53B  /  USS  Monterey 
transit  of  Bab  al  Mandab 


Figure  1 1 .  AN/SQS-53B  shallow  water  performance  at  beam  angles  of  010°  and  270°.  Data  collected  at  25  kts. 


Figure  12a.  AN/SQS-53B  shallow  water  performance  at  beam  angle  of  10°  showing  acoustic  returns  (bright 
areas)  from  the  shoal  areas  surrounding  Perim  Island. 


Figure  13.  Geo-rectified  AN/SQS-53B  sonar  image  forward  of  the  ship  as  it  heads  in  a  northward  direction. 
The  track  of  the  ship  is  traced  by  the  bright  or  whited-out  portion  of  the  imagery.  Note  that  when  the  ship 
changed  course,  the  look  direction  shifted  from  the  port  side  to  the  starboard  side  of  the  vessel  due  to  the 
receiver  beams  being  stabilized  (i.e.,  each  beam  remains  pointed  in  a  constant  direction,  relative  to  true  north, 
regardless  of  the  ship’s  course). 


Appendices 


1 .  Conversion  from  range  to  horizontal  distance  under  a  flat  bottom  assumption. 

2.  Specialized  AN/SQS-53B  Geo-rectification. 


1.  Conversion  from  range  to  horizontal  distance  under  a  flat 
bottom  assumption 


Let  Ti  be  the  range  for  sample  i  given  by  Vi  =  iAr  where  Ar  is  the  sample  spacing 
given  by 


where 


c  =  speed  of  sound 
At  =  samping  period 


For  a  set  of  samples  equally  spaced  in  range,  and  a  water  depth  z,  a  set  of 
corresponding  horizontal  distances,  Xj,  is  given  by 


However,  these  Xi  are  not  evenly  spaced  as  desired.  In  fact 


In  order  to  have  a  set  of  uniformly  spaced  horizontal  values  (spaced  at  Ar  to  agree 
with  the  long  range  case)  we  must  interpolate  between  the  Xi  to  get  a  set  x'^  which 
are  equally  spaced.  Let  Ax'  =  Ar  =  Choose  a  large  i  =  m  such  that 


1  - 


Ar 


<  e 


for  some  arbitrarily  small  e.  Values  of  the  acoustic  signal  at  the  new  sample 
points  are  obtained  through  interpolation  (including  the  special  case  of  "nearest 
neighbor”).  The  following  algorithm  can  be  used: 


^m  =  ^m 


where  m  is  deterinined  as  above. 


fori  =  - { 

A  =  ^i+i  - 


Find  j  such  that 
Then 


r'  = 

rj<r\<  Tj+i 


{^j+i  ~ 
(rj+i  -  r^) 


Air'iix';))  =  A{rj)  +  {r'i-rj) 

} 

where  A  (vj)  =  Aj  is  the  signal  amphtude  at  sample  j  at  range  Vj. 
1.  The  geometry  is  shown  in  Figure 


2.  Specialized  AN/SQS-53  Georectification 

Since  all  the  pings  for  a  given  beam  are  geostabihzed  and  thus  parallel  (for  fixed 
Sector  Center  direction),  simphfications  are  possible  in  performing  the  georectifi¬ 
cation  of  the  pings  The  basic  idea  is  to  construct  an  orthogonal  coordinate  system 
{x,y)  such  that  the  x-axis  is  parallel  to  the  ATth  beam  pings.  Assuming  that  the 
water  depth  correction  (flat-bottom  or  exact)  has  been  made,  the  indexing  of 
each  ping  is  modified  so  that  for  all  samples  (of  pings  •  •  •  j,j  -I- 1,  j  -h  2,  •  •  ^which 
have  the  same  x-coordinate  have  the  same  index,  i.  Figure  shows  the  concept. 
In  this  figure  the  {U,  V)  axes  are  geodetic,  e.g.  UTM  while  the  {X,  Y)  axes  are 
rotated  such  that  X  is  in  the  direction  of  the  Nth  beam.  Once  the  indices  of  each 
ping  is  adjusted,  we  construct  a  grid  with  a  uniform  spacing  along  the  y-axis  by 
interpolation.  The  resulting  set  of  pings  will  then  be  gridded  or  equally  spaced 
in  X  and  Y  although  the  spacing  in  X,  Ax,  and  in  Y,  Ay  will  not  necessarily  be 
the  same.  The  geographic  coordinates  of  the  four  corners  of  this  grid  are  known 
which  allows  for  a  trivial  georectification  using  an  internal  module  of  most  any 
image  processing  software  package.  The  geometry  is  shown  in  Figure 
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Figure  1.1: 
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2.1.  Along  ping  georectification 

Assume  the  ship  position  at  the  start  of  ping  j  is  given  by  {ujVj)  and  that  the 
angle  of  the  Nth  beam  is  ^  =  90°  -  Bn  Qs  shown  in  the  figure.  Let  the  pings 
be  indexed  {0, 1, 2,  •  •  • ,  j,  j  +  1,  •  •  • ,  </}•  Begin  the  procedure  by  translating  the 
{U,  V)  origin  to  (uo.’^^o)-  the  begiiming  of  the  first  ping. 

u'  =  u  —  Uq 
v'  =  V  —  Vo 

Then  rotate  the  axes  by  6  to  obtain 

u'  =  X  cos  6  —  y  sin  6 
v'  =  X  sin  9  +  y  cos  6 


cos  9  —  sin  ^ 

X 

i?' 

sin  9  cos  9 

.  y . 

The  new  coordinate  values  are  given  by 
or 


or 


X 

cos  9 

—  sin  9 

-1 

u' 

y 

sin  9 

cos  9 

v'  _ 

X 

y 


L 

J 

L 

cos  9 

sin^ 

■  u'  ' 

—  sin0 

cos  9 

v' 

or 


X  =  u'  cos  9  +  v'  sin  9 
y  =  —u'  sin  9  +  v'  cos  9 


or 


X  =  {u-  uo)  cos  9  +  {v  —  vq)  sin  9 
y  =  —  (u  —  Uq)  sin  9  +  {v  —  uq)  cos  9 
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In  the  new  (X,  Y)  system  we  find  the  ping-staxt  with  the  minimum  x-coordinate, 
j  =  M  and  the  maximum  x-coordinate,  j  =  Q  if  —90°  <  9  <  +90°.  (If  the 
pings  point  in  the  negative  x-direction  pick  the  maximum  and  minimum,  respec¬ 
tively).  We  then  translate  the  x-axis  to  Xm' 

X  =  x'  +  Xm 

y  =  y' 

With  these  transforms  the  locations  of  every  ping  sample  for  the  particular 
beam  N  should  be  positive.  Let  =  a(xij,yij))  be  the  value  of  the  ith  sample 
of  the  jth  ping  located  at  the  point  {xij,yij).  Let  aY  =  a(x'j,y'p  be  the  value  of 
the  sample  located  at  the  zth  sample  of  the  jth  ping  in  the  primed  coordinates 
(x',  y').  Recall  that  in  the  (x,  y)  coordinates  the  first  sample  occurs  at  the  location 
{xoj.yoj)  and  that  for  all  pings  except  for  j  =  M,  XQj  >  xqm-  Let  the  sampling 
interval  along  the  x-axis  be  Ax.  If  the  jth.  ping  begins  at  {xoj,yoj)  then  the 
coordinate  of  the  ith  sample  is  (xqj  +  iAx,yoj).  In  the  primed  system  let  the 
values  of  the  samples  for  x'^j  <  xqj  —  xm  be  zero.  Let 


^  Ax 


so  that  the  index,  kj,  corresponds  to  the  first  non-zero  value  of  ping  j. 

Thus  we  can  index  the  values  a[j  in  the  rotated  and  shifted  system  (x',  y')  as 
follows; 


^ij 


0  0  <i  <  kj 

^{i—kj)j  kj  i  ^  kj  +  Iji 


where  /max  is  the  mnnber  of  samples  in  the  original  (depth-corrected)  pings  (all 
equal).  If  we  require  that  each  ping  have  the  same  number  of  samples  then  we 
must  add  enough  zero  values  to  each  ping  to  give  a  total  of 


Tt  _  (^Q  ^m)  I  J 

-'max  ~  -'max 


With  this  notation  we  have 


0  0  <i  <  kj 

0-{i—kj)j  kj  i  kj  +  /max 
0  kj  +  /max  ^  —  -^max 
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2.2.  .Ping-toping  rectification 

The  gridding  along  the  y-axis  consists  of  choosing  a  Ay  and  calculating  a  new  set 
of  equally  spaced  pings  at 

ym  =  yo  +  mAy 

Since  the  are  already  aligned  it  only  remains  to  interpolate  in  the  y- 
direction.  The  algorithm  for  linear  interpolation  follows: 


for(m 

find  j  such  that  yj 
for  {i 


< 


0;  m  <  TUmax;  m  +  +{ 
yo  +  mAy  <  yj+i 

Oj  Z  ^  Zmaxi  ^  ”1  { 


^ij  ”1”  {V'm  yj) 


{o-'ij  Q^o+i)) 

{yj  ~  yj+i) 


} 

} 


where  the  zth  subscript  on  y  is  omitted  since  it  is  redundant. 


7 


